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O ( The colour fields created by the static tetraquark and pentaquark systems are computed in 

quenched SU(3) lattice QCD, with gauge invariant lattice operators, in a 24 3 x 48 lattice at 
i— i j3 =6.2. We generate our quenched configurations with GPUs, and detail the respective bench- 

manrks in different SU(N) groups. While at smaller distances the coulomb potential is expected 
to dominate, at larger distances it is expected that fundamental flux tubes, similar to the flux-tube 
^j- between a quark and an antiquark, emerge and confine the quarks. In order to minimize the po- 

tential the fundamental flux tubes should connect at 120o angles. We compute the square of the 



colour fields utilizing plaquettes, and locate the static sources with generalized Wilson loops and 
with APE smearing. The tetraquark system is well described by a double- Y-shaped flux-tube, 
with two Steiner points, but when quark- antiquark pairs are close enough the two junctions col- 
lapse and we have an X-shaped flux-tube, with one Steiner point. The pentaquark system is well 
described by a three- Y-shaped flux-tube where the three flux the junctions are Steiner points. 
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Figure 1: We compare (a) the tetraquark flux tube (or string) model, the elementary flux tubes meet 
in two Fermat points, at an angle of a = 120° to form a double- Y flux tube (except when this is 
impossible and the flux tube is X-shaped) with (b) Lagrangian density 3D plot for r\ = 8, ri = 14, 
presented in lattice spacing units (colour online). 

1. Motivation 

Multiquark exotic hadrons like the tetraquark and the pentaquark, different from the the or- 
dinary mesons and baryons, have been studied and searched for many years. The tetraquark 
was initially proposed by Jaffe [1] as a bound state formed by two quarks and two antiquarks. 
Presently several observed resonances are tetraquark candidates. The most recent tetraquark candi- 
dates have been reported by the Belle Collaboration in May, the charged bottomonium Z^(10610) 
and z£( 10650) [2]. However a better understanding of tetraquarks is necessary to confirm or dis- 
prove the X, Y and possibly also light resonances candidates as tetraquark states. 

On the theoretical side, the first efforts have been to search for bound states below the strong 
decay threshold [3, 4, 5, 6], as it is apparent that the absence of a potential barrier may produce a 
large decay width to any open channel. Recent investigations found that the presence of an angular 
momentum centrifugal barrier may increase the stability of the system [7, 8]. 

In the last years, the static tetraquark potential has been studied in Lattice QCD computa- 
tions [9, 10, 11]. The authors concluded that when the quark-quark are well separated from the 
antiquark-antiquark, the tetraquark potential is consistent with One Gluon Exchange Coulomb po- 
tentials plus a four-body confining potential, suggesting the formation of a double-Y flux tube, as in 
Fig. 1, composed of five linear fundamental flux tubes meeting in two Fermat points [6, 12, 13]. A 
Fermat, or Steiner, point is defined as a junction minimizing the total length of strings, where linear 
individual strings join at 120° angles. When a quark approaches an antiquark, the minimum poten- 
tial changes to a sum of two quark-antiquark potentials, which indicates a two meson state. This is 
consistent with the triple flip-flop potential, minimizing the length, with either tetraquark flux tubes 
or meson-meson flux tubes, of thin flux tubes connecting the different quarks or antiquarks [6, 8]. 
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Figure 2: Tetraquark Wilson loop as defined by Alexandrou et al , and by Okiharu et al 

Here we study the colour fields for the static tetraquark system with the aim of observing the 
tetraquark flux tubes suggested by these static potential computations. The study of the colour fields 
in a tetraquark is important to discriminate between different multi-quark Hamiltonian models. 
Unlike the colour fields of simpler few-body systems, say mesons, baryons and hybrids, [14, 15, 
16, 17], the tetraquark fields have not been previously studied in lattice QCD. 



2. Computing Fields with the Wilson loop and the Plaquette 

To impose a static tetraquark, we utilize the respective Wilson loop [9, 10] of Fig. 2, given by 
W 4Q = iTr(Mi^i 2 M 2 Li 2 ), where 

r>aa f ^ „abc ^a'b'c' nbb f ncc f 

K n — 2 1 2 ' 

r aa' 1 n abc n a'b ' c' r bb' r cc' si i\ 

The chromoelectric and chromomagnetic fields on the lattice are given by the Wilson loop and 
plaquette expectation values, 

<^(r)) = <P(r) M )-^f' r )y (2.2) 
\ i w/ \ v mi (w(n,r2,T)) 

(W)(n,r 2 ,T) \^ r W> 

where the jk indices of the plaquette complement the index i of the magnetic field, and where the 
plaquette at position r = (x,y,z) is computed at t = r/2, 

P^tv (r) = 1 - ^Re Tr [t/ M (r)U v (r + (r + v)Ul (r)] . (2.3) 

The energy (^f 7 ) and lagrangian (j£f ) densities are then computed from the fields, 

WO) = ^«E 2 (r)) + (B 2 (r)» , (2.4) 

(Jzf(r)) = i((E 2 (r))-(B 2 (r)» . (2.5) 
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Figure 3: Lagrangian density for Y2 = 14 and r\ from to 6. 
the Fermat points. The results are presented in lattice spacing 



The black dot points correspond to 
units (colour online). 



To produce the results presented in this work , we use 1121 quenched configurations in a 
24 3 x 48 lattice at /3 = 6.2. We also test that these configurations are already close to the continuum 
limit in a larger, 32 3 x 64 lattice. We present our results in lattice spacing units of a, with a = 
0.07261(85) fm or a~ l = 2718 ± 32 MeV. We generate our configurations in NVIDIA GPUs of 
the FERMI series (480, 580 and Tesla 2070) with a SU(3) CUDA code upgraded from our SU(2) 
combination of Cabibbo-Marinari pseudoheatbath and over-relaxation algorithm [18, 19]. 

To compute the static field expectation value, we plot the expectation value (£?(r)) or (#?(r)) 
as a function of the temporal extent T of the Wilson loop. In order to improve the signal to noise 
ratio of the Wilson loop, we use 50 iterations of APE Smearing with w = 0.2 (as in [16]) in the 
spatial directions and one iteration of hypercubic blocking (HYP) in the temporal direction. [20], 
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(d) Lagrangian Density 



Figure 4: Colour fields, energy density and Lagrangian density for r\ = 8 and r2 = 14. The black 
dot points correspond to the Fermat points. The results are presented in lattice spacing units (colour 
online). 



with ot\ = 0.75, &2 = 0.6 and = 0.3. At sufficiently large T, the groundstate corresponding to 
the studied quantum numbers dominates, and the expectation value of the fields tends to a hori- 
zontal plateau, for each point r determined by the plaquette position. For the distances r\ and r2 
considered, we find in the range of T G [3, 12] in lattice units, horizontal plateaux with a % 2 /dof 
E [0.3,2.0]. We finally compute the error bars of the fields with the jackknife method. 



3. The Tetraquark fields 

In our simulations, the quarks are fixed at (± r\ /2, — r2/2, 0) and the antiquarks at (± n /2, 7*2/2, 0) , 
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Figure 5: Wilson loop geometry, and preliminary result for the field, (Lagrangian density) for the 
static pentaquark. The results are presented in lattice spacing units (colour online). 

with r\ extending up to 8 lattice spacing units and Y2 extended up to 14 lattice spacing units, in order 
to include the relevant cases where Y2 > \fi>r\ . Notice that in the string picture, at the line r2 = V^n 
in our (n, r^) parameter space, the transition between the double- Y, or butterfly, tetraquark geom- 
etry in Fig. la to the meson-meson geometry should occur. The results are presented only for the 
xy plane since the quarks are in this plane and the results with z ^ are less interesting for this 
study. The flux tube fields can be seen in Fig. lb, 3 and 4. Theses figures exhibit clearly tetraquark 
double- Y, or butterfly, shaped flux tubes. The flux tubes have a finite width, and are not infinitely 
thin as in the string models inspiring the Fermat points and the triple flip-flop potential, but never- 
theless the junctions are close to the Fermat points, thus justifying the use of string models for the 
quark confinement in constituent quark models. 

We also compare the chromoelectric field for the tetraquark and the quark-antiquark system 
in the middle of the flux tube between the (di)quark and the (di)antiquark, and confirm that the 
tetraquark flux tube is composed of a set of fundamental flux tubes with Fermat junctions. 

4. Foreword 

• The flux tubes remain interesting in Lattice QCD, our results support the string model of 
confinement, in particular for the tetraquark static potential 

• The mixing between the tetraquark and meson-meson flux-tubes is small, which may con- 
tribute for narrower tetraquark resonances. 

• We are now studying in more detail other flux tubes, as for the pentaquark in Fig. 5. 
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